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Abstract. - Within an Onsager-like density functional theory we explore the thermodynamic 
and structural properties of an isotropic and nematic fluid of hard needle-like colloids in contact 
with a hard substrate coated with a soft short-ranged attractive or repulsive layer. As a function 
of the range and the strength of the soft interactions we find wetting and drying transitions, 
a pre-drying line, and a symmetry-breaking transition from uniaxial to biaxial in the wetting 
and drying film. 



Whereas bulk liquid crystal phases of suspensions of colloidal hard rods have been essen- 
tially understood due to Onsager's work in the 1940's [1,2] and simulations [3] and density 
functional theories [4-6] in the late 1980's, their surface and interfacial properties are subject of 
ongoing study, not only experimentally but also theoretically and numerically. Good progress 
was made during the past decade in the theoretical study of planar free isotropic-nematic (IN) 
interfaces, e.g. it is known by now that the nematic director in the thermodynamically stable 
state is parallel to the interfacial plane [7-9], that complete wetting of the IN interface by 
another nematic phase occurs near the triple point in binary rod mixtures [10], and that the 
one-particle distributions in the IN interface are only weakly biaxial [11-13]. The effect of ex- 
ternal substrates on suspensions of hard rods has also been studied. For a planar hard wall, for 
instance, evidence of complete wetting of the wall-isotropic (WI) interface by an intervening 
nematic film was provided by theory [14] and simulations [15]. Other studies were concerned 
with properties of a hard-rod fluid in contact with a "penetrable" wall, which restricts only 
the translational degrees of freedom of the rods [16]. It was shown that such a wall favors 
homeotropic anchoring of the nematic director, and that the WI interface exhibits complete 
wetting by the homeotropically aligned nematic phase upon approach of IN coexistence [8,17]. 
The common feature of all these studies is that the chosen wall potential does not allow to 
control the degree of surface nematic order. This is in contrast with Landau-de Gennes the- 
ory, which predict rich surface phase diagrams for liquid crystals [18-20]. A drawback of this 
formalism is, however, that the effects of particular surface-particle interactions are hidden in 
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expansion coefficients that are not always easily traced back to the microscopic details [20,21]. 
Therefore, it is interesting and relevant to construct surface phase diagrams, for a "simple" 
fluid of monodispcrse hard rods, in terms of parameters of the surface potential. 

Calculating surface phase diagrams of colloidal liquid crystals from a microscopic theory 
is technically involved and computationally expensive, as density profiles p(z, 6, ip) are to be 
determined for a given wall potential V{z, 9, (f>), with z the distance from the wall and and ip 
the polar and azimuthal angle of the orientation of the rod. Building on earlier studies [22,23], 
only recently a study appeared of the surface properties of a hard-spherocylinder fluid in 
contact with a model substrate, composed of a "penetrable" wall and a short-ranged repulsive 
or attractive tail [17]. Several wetting transitions and a transition from homeotropic to planar 
anchoring were observed at different strengths of the wall potential. This can be understood 
qualitatively, since the "penetrable" wall with the attractive tail favors homeotropic alignment 
of the nematic director, and for the strong repulsive tail it resembles the hard wall such that 
the stable director changes to planar. However, the complexity of the model prohibits to 
explore the full parameter space. 

In order to be able to map a complete surface phase diagram we sacrifice the full rotational 
symmetry of the rods and consider the simpler Zwanzig model, with a restricted number 
of allowed orientations [24]. This model exhibits a strong first-order IN transition in the 
bulk [24] , and the orientation of the nematic director parallel to the IN interface was found 
to be thermodynamically favorable [25], similar to continuous rods. In contact with a hard 
wall the WI interface is completely wet by the N phase upon approach of the IN coexistence, 
whereas the WN interface remains partially wet [25], i.e. in such a geometry the model also 
shows behavior similar to continuous rods. Although some limitations of the Zwanzig model 
are known [26,27], we adopt it for the present study for reasons of its numerical simplicity. 
We explore the complete surface phase behavior of a Zwanzig hard-rod fluid in contact with 
a substrate which consists of a hard wall and a short-ranged attractive or repulsive tail. In 
contrast with the "penetrable" wall it puts the infinite potential barrier at the origin regardless 
the orientation of the rod. 

We consider an inhomogeneous fluid of rectangular hard rods of length L and diameter D 
(L 3> D) in a macroscopic volume V at inverse temperature f3 and chemical potential /i. The 
rod orientations are restricted to the three mutually perpendicular directions hi representing 
x, y, z for i = 1, 2, 3, respectively. The position r of the center of mass of a rod is continuous. 
The grand potential functional of the fluid in an external potential Vi (r) with the one-particle 
distribution functions Pi(r) can be written, within the second virial approximation, as [25] 

0CI\P]=Y, I dv Pi {v)(Hp i {v)v]-l-0 l ji + (3V i {v))-\ [ ^r'.^(r;r')ft(r)^(r'),(l) 

»=i »,j=i 

with v = L 2 D, and fij(r; r') = —1, (0) the Mayer function of two (non-)overlapping hard rods 
with orientations hi and hj and center-of-mass coordinates r and r'. Qualitatively, the second 
virial functional predicts the bulk IN transition and the adsorption and wetting properties 
near a planar hard wall in agreement with simulations of freely rotating rods [25]. 

The external potential Vi(r) consists of a single planar hard wall in the plane z = (normal 
z) and a generic soft "tail" oc exp(— kz) acting on each rod segment, and is written as 

{PAexp[— kz] 2>0,i = l,2; 
f3AFexp[-Kz] z>±,i = 3; (2) 
oo otherwise, 
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where the form factor F = sinh(KL/2)/(KL/2) > 1 arises from integration over the contour 
of the rod. For A > such a potential could , under certain circumstances, be caused 
by a polymer coating, with A a measure for the planar density and for the radius of 
gyration [28]. In the present study we avoid discussions on the microscopic origin of Vi(r), 
and explore the generic properties of the surface phase diagram as a function of the "contact" 
potential A and the decay constant k. Our primary interest is in the parameter regime 
j3A < 10, since larger values of (3A correspond essentially to a shift of the hard wall along the 
z-axis. 

Throughout we assume translational invariance in the xy plane, such that the dimensionless 
densities Ci(z) — Pi(v)v only depend on the distance from the substrate z. The Euler-Lagrange 
equations SD,[{p}]/Spi(r) = can then be written as [25] 

0fx = ln Cl (z) + 0Vi(z) + 2c 2 (z) + 2c 3 (z), 

Pfj, = lnc 2 (z) + (3V 2 (z) + 2ci(» + 2 c 3 (z), (3) 
j3p = lnc 3 (z) + 0V 3 (z) + 2c x {z) + 2c 2 (z), 

with the averaged densities Ci(z) — (l/L) J_ L , 2 dz'ci(z + z'). We solve Eqs. iteratively 
for a given A and k at a fixed chemical potential p. In all numerical calculations we use an 
equidistant z-grid of 40 points per L. Convergence is assumed when the relative difference 
between the results of the subsequent iterations is smaller then 10 -10 for all values of z in the 
grid. Such accuracy is required to avoid dependencies of the results of the calculations on the 
initial guesses. Additional checks, performed with 80 points per L, yield virtually identical 
results. Once determined, the equilibrium profiles can be inserted into the functional to obtain 
the equilibrium value of the grand potential 

^ = ^E j dzc i {z)(Hc i {z)]-2-f3ix + (iV i {zj), (4) 

Note that f2 = — pV for a bulk system in a volume V, with p = p(p, T) the pressure. In 
the presence of a planar surface or interface of area A we have f2 = — pV + jA with 7 = 
7(/i,T) the surface or interface tension. The results of the calculations can be conveniently 
represented in terms of the total density c(z) = ^2 i= i Ci(z), the nematic order parameter 
s(z) — (03(2) — i(ci(z) + 02(2))) / c(z) , and the biaxiality A(z) = (ci(z) — C2(z))/c(z). The 
(dimensionless) adsorption near a planar substrate is defined as r(/x) = dz(c(z) — ct,)/L, 
where q, = Cb{p) is the bulk density far from the substrate. In Fig. ^a,b,c) we show p — A 
surface phase diagrams for several ranges of the wall potential kL = 6 (a), nL — 2 (b), and 
kL = 1 (c). The phase diagrams are very rich, therefore we discuss separately their parts 
which correspond to the low (/i < /i/at) and high (// > pin) values of the bulk chemical 
potential. 

For p < pin all three phase diagrams of Fig. ^a,b,c) exhibit a uniaxial-to-biaxial (UB) 
surface transition at fi = pub(A) < pin for small enough A < Ay B . This implies that 
only sufficiently repulsive walls do not show the UB transition, such that A(z) = even 
upon approach of IN coexistence. The critical value satisfies {iA\j B > 0.5 in all cases, and 
increases with kL (Fig.QJd)). The UB transition is found to be continuous for all investigated 
combinations of kL and (3 A. This is consistent with the findings of Ref. [25], where the case 
(3 A = (hard-wall) was considered. 

For all studied kL we only recover the hard- wall complete wetting phenomenon [25]) for 
sufficiently weak repulsions A < A w , whereas strong repulsions A > A w give rise to partial 
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Fig. 1 - (a,b,c) Surface phase diagrams of the Zwanzig hard-rod fluid in contact with a soft wall in 
terms of the contact potential j3A and chemical potential for different range of the wall potential 
kL — 6 (a), kL = 2 (b), and kL = 1 (c). In all three cases we distinguish a uniaxial (U) and biaxial 
(B) phase for /i < [iin, and a wetting transition (*) at /i = fiiM that separates complete wetting of 
the wall by a nematic film at A < A w from partial wetting A > A w . In (b) and (c) we also see a 
predrying line at \i > [im, and a drying transition at /3fj, = PfiiN, that separates complete drying by 
an isotropic film at A > Ad from partial drying A < Ad- The inset in (a) illustrates the logarithmic 
divergence of the adsorption T with undersaturation e = 1 —Cb /cj for A < A w , and that in (c) shows 
several biaxiality profiles A(z) for A = 0.55 > Afj B for several values of the bulk chemical potential 
£ti,2,3 > fiiN- (d) Generic surface phase diagram for a fluid of Zwanzig hard rods at IN coexistence 
in contact with the "soft" wall. Solid lines correspond to the wetting (□), drying (o), and the UB 
(o) transitions; the end point on the line of drying transitions is indicated by (*). The dashed lines 
indicate k£i,n ~ 1, with the bulk correlation lengths £i/L — 0.54 and £n/L = 0.35 of the bulk I and 
N phase, respectively. 



wetting. This is illustrated in the inset of Fig. 1(a), where the adsorption T diverges logarith- 
mically with the undersaturation e = 1 — for A < A w ~ 1.1, and remains finite 
for A > A w . Here cj is the coexisting isotropic bulk density. This implies that a wetting 
transition takes place at A = A w . Interestingly, the numerical value A w ~ 10 -2 is seen 
to be very small for longer-ranged potentials, kL < 2, i.e. extremely weak but sufficiently 
long-ranged repulsions strongly reduce the complete- wetting regime. Shorter-ranged substrate 
potentials exhibit a larger complete wetting regime, e.g. A w =~ 1.1 when kL — 6 (see also 
Fig^d)). The existence of a wetting transition at A = A w follows also from comparisons of 
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Fig. 2 - (a) The relative surface tension difference R(A) at the WI and WN interfaces for a wall 
potential with a decay given by kL = 6. The vertical dashed line /3A m — 1.15 indicates the repulsion 
amplitude at which the wetting transition takes place. The inset shows the excess adsorption T(A) at 
coexistence, (b) Excess desorption F(e') at the interface between the N phase and the soft wall (kL = 
2) as a function of the oversaturation e' = cs/cjv — 1 for several different values of flA — 0.721(1 + S). 
The vertical dashed lines connect coexisting points at the predrying line. 



the corresponding interfacial energies. For convenience, we define the relative surface tension 
difference 

= 1wi (A)-1wn{A)^ 

JIN 

with jwi(A) and Jwn(A) the tensions of the WI and WN interface, respectively, both at 
A« = fJ-iN, and (3jinLD — (2.8027 ± 0.0001) x 10~ 2 the surface tension of the free planar IN 
interface [25]. In Fig. Efa) we show R(A) for kL = 6. It can be seen that R(A) = 1 (with 
relative accuracy of 10 -3 ) if f3A < f3A w = 1.15, and R{A) < 1 if A > A w . This is a clear 
thermodynamic evidence for a transition from partial to complete wetting at A = A w . 

The study of the order of the wetting transition represents significant numerical difficulties. 
Our results suggest that it is likely a first-order wetting transition, because (i) R(A) appears 
to approach R(A W ) = 1 with a finite slope when A j. A w , (ii) metastable nematic films of finite 
thickness can easily be generated at n = (iin for A < A w , and (iii) the adsorption T(/iin) at 
saturation appears to show a discontinuous jump from a finite value for A > A w to an infinite 
value for A < A w (see insets Fig.^a) and^a)). One expects, for reasons of continuity, that 
a first-order wetting transition is accompanied by a first-order prewetting transition, from a 
finite thin to a finite thick film off coexistence. Unfortunately, we have not been able to find 
this prewetting transition, despite considerable efforts. We speculate that its critical point is 
too close to the wetting transition to be detected. 

We now consider the /it > fiiN part of the phase diagrams, where the N bulk phase has 
director x, i.e. planar. We have checked that a homeotropic orientation of the director is 
metastable in all the cases studied. Intuitively, one expects a low-density isotropic film close 
to a highly repulsive substrate {(3 A 3> 1), which may but need not become of macroscopic 
thickness when fi J, /Ujjv- For (3 A = it is known that T < but finite [25], a situation that we 
refer to as partial drying. Our numerical results for kL = 6 show that T remains finite when 
A 1 I MiJV for any (3 A < 10. For (3A > 10 the density of rods close to the hard wall (z = 0) 
is so small that its position can be shifted to some z = z > without affecting the profiles 
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Ci(z) significantly, while the effective amplitude is reduced from A to Aexp(— hzq) < 10. We 
can conclude that the soft potential with kL = 6 is so short-ranged that the partial drying is 
obtained for any value of (3 A. Fig. 1(a) is therefore featureless for [i > ii/Af- 

The situation is more interesting for the longer-ranged cases kL = 2, 1, where we do 
find a drying transition at (3Ad — 0.72,0.31 in Fig. ^b), (c), respectively. When /i j. /xjjv 
partial drying is observed for A < Ad, and complete drying by an isotropic film for A > Ad- In 
Fig.|^b) we show the excess adsorption r(e') as a function of oversaturation e' = c^/cm — 1 for 
various amplitudes of the wall potential with kL = 2. At low values of the contact potential 
(A < Ad) the desorption remains almost constant upon approach of the IN coexistence. 
Calculations show that R(A) > —1 in this regime, corresponding to partial drying. The 
behavior changes qualitatively for 0.72 < [3A < 0.95, where r(e') exhibits a discontinuity at 
some e' > 0, which we associate with predrying. For (3 A ~ 0.95 the jump of the desorption 
takes place at e' = 0.05, and its magnitude is (vanishingly) small, i.e. this is the critical 
predrying point indicated by (*) in Fig. ^b) and (c). The jump of the desorption at the 
predrying line increases upon decreasing (3 A towards (3 Ad, and its location shifts to smaller 
e', until it diverges at the drying transition at (3A — f3Ad at e' — > 0. For (3A > 1.0, i.e. 
above the critical predrying amplitude, the desorption grows continuously upon approach of 
coexistence. Studies of the surface tensions at coexistence show that R(A) = —1.000 ± 0.001 
for all A > Ad, consistent with a complete drying of the WN interface by an intervening 
isotropic film. However, the profiles (not shown) do not exhibit a clear bulk film even at 
e' = 10 -5 , which implies that the asymptotic limit e' — ► has not been reached yet such that 
the regime of logarithmic growth of — r with e' cannot properly be identified in Fig. Gib). 
Interestingly, for kL < 1.5 we find that Ad < Ay B , such that for A > Ad the (thick) isotropic 
film that grows in between the wall and the nematic bulk phase can exhibit a UB transition. 
In other words, the UB transition at fi < ///at extends continuosly to the thick isotropic films 
at n > hin, as illustrated for kL = 1 in Fig. E^c). Here the line of UB transitions terminates 
at (3Ay B — 0.63, where it crosses the predrying line. We note, finally, that no new physics is 
to be expected for larger A in FigQIb) and (c), since the regime of approximate scaling of the 
density profiles (associated with a shift of the hard wall and a renormalization of the contact 
potential) at large A is observed already for (3 A > 1.5. 

All our findings are summarized in Fig. HTd). which represents the generic surface phase 
diagram in the (3A — kL coordinates for a fluid of Zwanzig hard rods at IN coexistence 
(/i = (J,in) in contact with the "soft" wall. We distinguish lines of the wetting (□) and drying 
(o) transitions. The UB transition (o) at the WI interface occurs for all kL. It also takes 
place at the WN interface in the thick isotropic wetting film, i.e. below the intersection 
of the lines of the drying and the UB transitions. Note that the growth of the wetting 
film is determined by the correlation length £j (£jv) of the coexisting bulk isotropic and 
nematic phase, and no longer by when kT x < £j (&v)- However, this crossover does 
not introduce any additional structure to the surface phase diagram. For large kL the UB 
line shows an exponential divergence (3 A ~ exp[/cL], associated with the renormalization of 
the contact potential in the scaling regime (3 A = (3Aq exp[K,Lzo/L]. The line of the wetting 
transitions is seen to be linear with kL > 4. This is associated with the fact that for the 
steep potentials density variations within the decay length of the potential are rather small, 
and can be neglected. Therefore, its contribution to the grand potential is proportional to 
((3A/ L)c(0) J dz exp[—Kz] — ((3A/ kL)c(0). Since the wetting transition for large kL occurs at 
similar values of the surface free energy, {(3A W / 'kL)c(O) = const, and hence f3A w ~ kL. 

We conclude that the surface phase behaviour of hard rods with soft walls is extremely 
rich, at least on the basis of the present analysis within the Zwanzig model. We find that 
complete wetting of the WI interface by a nematic film, as found for a planar hard wall, 
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is very sensitive to weak but long-ranged wall repulsions, and that strongly repulsive walls 
show complete drying of the WN interface by an intervening isotropic film if the range of 
the potential is long enough. It is of interest to extend the present study to other related 
systems, e.g. freely rotating rods or platelets. Another related interesting question concerns 
the required nature of the surface potential that can lead to a stable anchoring transition in 
a colloidal-rod fluid [23]. Work along these lines is in progress, and hopefully stimulates the 
study of colloidal analogues of liquid crystals near interfaces. 
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